Source of Acquisition 
NASA Goddard Space Flight Center 


Predicting the Starquakes in PSR J0537— 6910 

J. Middleditch 1 , F. E. Marshall 2 , Q. D, Wang 3 , E. V. Gotthelf 4 , and W. Zhang 2 

jon@lanl.gov 

ABSTRACT 

We report the results of more than seven years of monitoring of 
PSR J0537— 6910, the 16 ms pulsar in the Large Magellanic Cloud, using data 
acquired with the Rossi X-ray Timing Explorer. During this campaign the pul- 
sar experienced 22 sudden increases in frequency (“glitches” - 21 with increases 
of at least eight /iHz) amounting to a total gain of over six parts per million of 
rotation frequency superposed on its gradual spindown of v = —2 x 10 -1 ° Hz s _1 . 
The time interval from one glitch to the next obeys a strong linear correlation 
to the amplitude of the first glitch, with a mean, slope of about 400 days per 
part per million (6.5 days per /zHz), such that these intervals can be predicted 
to within a few daysman accuracy which has never before been seen in any other 
pulsar. There appears to be an upper limit of ~40 /zHz for the size of glitches 
in all pulsars, with the 1999 April glitch of PSR J0537— 6910 as the largest so 
far. The change of its spindown across the glitches, Az>, appears to have the 
same hard lower limit of -1.5 x 10~ 13 Hz s -1 , as, again, that observed in all other 
pulsars. The spindown continues to increase in the long term, v = -10~ 21 Hz s -2 , 
and thus the timing age of PSR J0537— 6910 (— 0.5z'z>~ 1 ) continues to decrease 
at a rate of nearly one year every year, consistent with movement of its magnetic 
moment away from its rotational axis by one radian every 10,000 years, or about 
one meter per year. PSR J0537-6910 was likely to have been born as a nearly- 
aligned rotator spinning at 75-80 Hz, with a \i>\ considerably smaller than its 
current value of 2xlO -10 Hz s -1 . Its pulse profile consists of a single pulse which 
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is found to be flat at its peak for at least 0.02 cycles. Glitch activity may grow 
exponentially with a timescale of 170 years viz {{vv)c Ta b)~ l in nil young pulsars. 

Subject headings: pulsars : neutron — stars:individual (PSR J0537— 6910) — X-rays:stars 


1. Introduction 

The X-ray pulsar, PSR J0537— 6910 (hereafter J0537), in the 30 Doradus star formation 
region of the Large Magellanic Cloud (LMC), was discovered serendipitously in an observa- 
tion of Supernova 1987A with the Rossi X-Ray Timing Explorer (RXTE) in a search for its 
pulsar remnant (Marshall et al. 1998; Middleditch et al. 2000). Associated with the 4,000- 
year old LMC supernova remnant, N157B (Wang & Gotthelf 1998), with a rotation period 
of only 16.1 ms (62 Hz), it is the fastest rotating young pulsar known, the next fastest two 
young pulsars being the Crab pulsar, with a period of 33 ms (30 Hz), and PSR B0540-69, 
with a period of 50 ms (20 Hz). With a spindown rate near —2 x 10 -10 Hz s -1 , J0537 has one 
of the highest energy loss rates known, of ~ 5 x 10 38 /4 5 ergs s -1 , where J45 is the moment of 
inertia of the neutron star (NS) in units of 10 45 g cm 2 , but it has so far only been detected 
in the X-ray band despite several radio and optical searches (Crawford et al. 1998; Mignani 
et al. 2000). 

Pulsars are known for their very stable rotation and small spindown rates, and their 
spindown histories can be used to constrain models of their emission mechanisms. Charac- 
teristic pulsar ages, r c , may be estimated, assuming that their spindowns are dominated by 
magnetic dipole radiation, v a u 3 , as r c = — 0.5pr> _1 , where v is the rotation rate in cycles 
s -1 (Hz), and v is its (almost always negative) time derivative, or spindown, in Hz s -1 . J0537 
has a young characteristic timing age of 5 x 10 3 years, consistent with the 4,000-yr estimated 
age of N157B, though we have showed this quantity to be actually decreasing with time 
(Marshall et al. 2004, hereafter M04). 

Many pulsars have also been shown to undergo sudden discontinuities (usually increases) 
in their rotation rates known as "glitches,” (see, e.g., Alpar et al. 1993,1996; Hobbs et 
al. 2002; Johnston et al. 1995; Jones 2002; Lyne, Pritchard, & Smith 1993; Lyne et al. 2000; 
Pines & Alpar 1985; Ruderman, Zhu, k Chen 1998; Shemar & Lyne 1996; Wang et al. 2000; 
Wang et al. 2001a; and Wong, Backer, & Lyne 2001). Although glitches have been detected 
in both young and old pulsars, they are predominately found in the younger ones (r c < 10 5 
yr, Urama and Okeke 1999, Lyne, Shemar, & Smith 2000). Glitches are thought to occur 
when angular momentum is transferred from a more rapidly rotating component of the NS 
to the outer crust (see also Anderson & Itoh 1975; Franco, Link, & Epstein 2000; and §5.1 
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for alternative models). 

The increase in rotation rate shows up within a few minutes because the NS star mag- 
netic field is thought to be fixed in the crust. Thus, in addition to the continuous spindown 
of a pulsar, glitches can be used to reveal details about the NS equation of state and internal 
structure, which would otherwise remain hidden (see, e.g., Datta & Alpar 1993; Link, Ep- 
stein, & Lattimer 1999). If, for example, erratic timing activity was found to always precede 
the large glitches, as seems likely from this work, it would help to eliminate mechanisms 
which rely on truly sudden onsets as plausible causes for glitching, such as sudden events 
involving the unpinning of vast numbers of superfluid vortices within the crust, triggered by 
mechanisms other than cracking in the solid crust. 

Below we report (§2) on more than seven years of observations of J0537, including more 
than four years of new observations, the whole of which argues that the time interval from 
one glitch (the first) to the next glitch (the second) in J0537 is strongly correlated to the 
amplitude of the first glitch (§3), a pattern that is similar to that of large quakes within the 
crust of our own planet.. In §5 we compare the behavior of J0537 with that of the population 
of glitching pulsars and discuss how its glitch/time correlation helps to discriminate between 
different proposed mechanisms triggering their glitches. Section 6 lists our conclusions. 


2. Observations 

The data were obtained with the Proportional Counter Array (PCA) on board the 
RXTE observatory, as described in detail in M04. The PCA is sensitive to X-rays in the 
2-60 keV band and has moderate spectral resolution (A E E~ l ~ 18%). Each event is 
time-tagged on the spacecraft with an uncertainty of less than 50 ps (Rots et al. 1998). The 
background consists of unrejected charged particles and X-rays from LMC X-l and other 
cosmic X-ray sources. Only events in the first xenon layer in channels 5-50 (~3-20 keV) 
were included in this study. A log of RTXE observations of J0537 is given in Table 1. 
The observations are separated by 22 glitches into two incomplete and 21 complete groups. 
The number and epoch range for observations for which we have a timing solution (i.e., 
an unambiguous cycle count linking its pulse time of arrival [TOA] to those of at least one 
neighbor) are listed under “Phased,” and the number and epoch range for the subset of. these 
which were used to derive an ephemeris are listed under “Fit.” Our monitoring system of 
J0537 continued as described in M04 except that, starting in late 2003 additional closely 
spaced observations were scheduled during the few week time period during which the next 
glitch was anticipated. 
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3. Data Analysis and Results 
3.1. The Pulse Profile 

We first corrected the photon arrival events to the Solar System Barycenter using the 
source position (J2000) ofa = 05 /ir 37 m 47,36 s , 6 = -69°10'20.4" (Wang et al. 2001b). Times 
are given in barycentric dynamical time (TDB). Other details are as given in M04. For 
each observation we determined the pulsar’s frequency and phase at the time of the first 
event (a nominal spindown, z>, of -1.988xlO -10 Hz s -1 was used throughout these processes). 
The pulse shape was fitted at first to an ad hoc downward parabola with a half width at 
zero height of 0.072 cycles, but later to a master pulsar profile (MPP) derived from fitting 
folded data from many observations to a Breit-Wigner function (see, e.g., Nelson et al. 1970). 
The uncertainty in phase was determined from Monte Carlo simulations to be ~0.04 SNR -1 - 
cycles, where SNR is the signal-to-noise ratio of the single peak of the folded J0537 pulse 
profile, including noise fluctuations due to background and the pulsar signal itself. For the 
parameters of our observations (listed in M04) only about 0.5% of the count rate is due to 
the pulsar. However, these events account for six times as much variance per event than 
does the background, due to their concentration within the pulse profile shape. 

Although several more iterations of this fitting made the next MPP sharper, eventually 
this process actually resulted in lower SNRs, at least until we utilized the frequencies as 
derived from the timing solutions (see §3.2 below), which were much, much more accurate 
than could ever be achieved by fitting a single observation, to derive improved phases from 
a simple, one dimensional Newton-Raphson fitting process. The resulting MPP appeared 
to be wider at the top of the peak than any previously derived MPPs, and could no longer 
be well fit by any Breit-Wigner (BW) formula. Accordingly, we subtracted a Gaussian of 
a given center, narrow width, and small amplitude, from the BW formula, and this fit our 
first timing solution MPP well. Another MPP was generated by fitting to this modified BW 
(MBW), and this appeared to be flat at the top for at least 0.02 cycles. This master pulse 
profile and the MBW used to produce it are plotted in Figure 1. 

Further iterations, with MPPs generated using the MBWs made by fitting to the MPPs 
from the previous iteration, could still be well fit by an MBW, but also produced MPPs that 
were progressively more “two-horned.” These MBWs also produced progressively worse 
SNRs when fit to the individual runs, a possible side effect of the interaction of statistics for 
the low SNR individual observations from which the MPP is produced, with the maximiza- 
tion process for the fit response. The MPP drawn in Fig. 1 is given by: 

L(<f>) = 1 ~ ~ -- y -- ? 'e~ C502 (l - c 6 e -C7 ^ - ^ 2 ), 

1 - (p{C 3 - (pCi) 


( 1 ) 
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where —0.5 < <p < 0-5- The three parameters that characterize individual runs were deter- 
mined by fitting a function of the form aL(<p — <j> 0 ) 4- b, where a gives the signal amplitude, 
b the background, and (p 0 defines the phase of the first pulse arrival in the observation time 
interval. A time of arrival (TOA) was obtained self-consistently by propagating the time 
of the first pulse peak arrival (derived from the measured pulse phase and the observation 
start time) , by an integral number of cycles to some epoch near the center-of-gravity of the 
observation interval, using the nominal i>, together with the very accurate frequency given by 
the timing solution, or lacking that, the measured frequency for the individual observation. 
The constants, Ci- 7 , ip, a, b, and <f>o are given to full accuracy in Table 2. The phase, <p 0 , like 
a and b, is in this case, for the MPP of Fig. 1, rather than any individual observation, and 
together with ip, ci_ 7 , a, and b, makes a total of eleven parameters. 


3.2. The Timing Solutions 

The method of determining the timing solutions was the same as has been frequently 
employed for pulsars and other periodic phenomena. Per M04, the initial time gap between 
two observations (~6-12 hours) was only a few times longer than the total time spanned by 
each observation, (5,000-20,000 s), usually consisting of 2-4 contiguously observed segments 
and their intervening gaps. This guaranteed that the average of the individual frequencies 
measured for two observations could reliably determine the mean frequency between the 
two. The time gap between the successive observations was then increased geometrically. 
A value for z>, close to those already determined from fits to other segments, or, lacking 
that, the long term frequency history of J0537, was used to count cycles over the longer 
time gaps. If the estimate of the number of cycles between TOA’s from two consecutive 
observations was within a few hundredths of an integer, and the measured frequencies of the 
two observations were consistent with the same frequency and the small correction for i>, one 
could be confident that the cycle-counting was unambiguous. 

Continued successes at cycle counting in extending the baseline of the timing solution 
quickly reduced the uncertainty of its validity to a vanishingly small value. When the time 
baseline of the observations and the gaps between them was sufficiently long, v could be 
reliably determined by the fitting process. Eventually, v could be determined when the time 
baseline exceeded 100 days. The interval at which the gap size was reset to ~6-12 hours 
was 60 days, and the whole timing procedure could be started over, if necessary. The cycle- 
counting process was extended over a longer and longer time baseline, until a glitch made 
further extension of the unambiguous cycle-counting impossible. 

A sum-squared minimizing routine was used to fit the data segments to a power law 
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model, involving the pulsar braking index, n, in v cc —v n . As an example, the frequency 
predicted by the power law relation at time t, can be, for the purposes of discussion, more 
conveniently expressed as: 

v{t) = v(t 0 )(l - (*>(io)/ v{to))(n - 1 )(t - to))~ 1/(n-1) , (2) 

where to is the epoch of the fit, and the quantity, — z>(t 0 )^(to)~\ is usually more concisely, 
though not as conveniently for the purposes of discussion, expressed as 7(^0) • As a check 
on the power law parameters and their errors, v, u, and i>, were derived, using simple linear 
algebra, from linear polynomial fits to the TOA’s and cycle numbers of points belonging to 
the individual data segments: 

mi - 4>o + v(toa,i — to) + v{toai — t 0 ) 2 /2 + i>(toa,i - 1 0 ) 3 / 6, (3) 

where m* is the integral number of cycles from time, to - to tooj, <j > 0 is the pulse phase 
at time, to, and the need for <p 0 can be eliminated by a proper choice of to. For another 
check, the power law fits were re-initialized with the equivalent parameters derived from the 
polynomial fits. The difference between the 3rd order polynomial fits using u, u, and i>, and 
the power law fit that uses the equivalent parameters, consisting of terms of higher time 
derivatives than u, was insignificant. 

In practice, when the data segments only contain a few points, they are first fit using 
only a fiducial phase or TOA, and v, with and n fixed, then with only n fixed, and 

then with all four parameters. The parameters of the fits and their errors are listed in Table 
3, along with an epoch of pulse peak arrival. In four cases of short data segments (5, 8, 
13, and 17), and three segments with >100% errors on n (1, 6, and 11), n was fixed at 3.5 
without much increase in chi-square. In three more segments (10, 14, and 20) with large 
errors on n and large chi-squared increases when n was set to 3.5, n was fixed at values hear 
to what the four-parameter fit would have produced (31.5, 35.0 and 19.0) . Of this last three, 
segment 14 showed evidence that the last point was high (see 3.4 and Figure 4). Segment 20 
may have had both a small, but long-lasting recovery (the last and highest point of Segment 
19 precedes it by only three days), and a later pre-glitch creep, as its favored braking index 
near 19 can only be diminished when observations were removed from near its beginning and 
its end, though more observation points would have been needed to be certain. This situation 
is likely to become more common as further segments with more near-glitch observations are 
obtained. 
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3.3. The Frequencies, and Spindowns 

The errors on the frequency measurements of the individual observations were also 
calibrated by Monte Carlo simulations, and found to be [0.0418,0.045 cycles] (SNR*T r7n5 ) _1 , 
for the upper and lower la frequency errors, respectively. Here, T rms is the square root of 
the mean squared moment in time for all of the events. For a continuous (i.e., unbroken) 
observation of duration T, T rms is (2V3)' 1 T, as long as the signal and background count 
rates don’t change with time. A simple way to understand these errors is that the phase error 
is the half width at half maximum, 0.044 cycles (Fig. 1) , divided by the SNR. The frequency 
error is then that error divided by T rms (see, e.g., Ransom, Eikenberry, &; Middleditch 
2002). Figure 2 shows the individually measured frequencies of J0537 with a convenient v 
subtracted, along with the lines/curves representing the frequency trends from the timing 
solutions. 

The 22 glitches, ranging in size from 0.016 to 0.68 parts per million (ppm), are labeled 
in Fig. 2. The mean rate of decrease, —1.9760 x 10 -10 Hz s -1 , is shown as the oblique dashed 
line. All .of the comments on the behavior of the glitching made in -M04 apply as well to 
this new total data set, now a factor of three times larger. The first glitch is still the largest 
of the 22, and a measurement of the braking index, n, characteristic of the actual physical 
mechanism which slows the rotation of J0537, appears to be as elusive as ever. However, the 
braking index of J0537 softens from 10.8 and 6.6 for the two halves of the long data segment 
2, and there is some evidence for the same softening in most other long data intervals. For 
the moment though, it is impossible to know if this is inconsistent with the behavior of the 
high braking indices in non-glitching pulsars found by Johnston and Galloway (1999). 

The history of z> is shown in Figure 3. The points are derived from timing solutions 
of subsegments of data with at least three observations, spanning at least 1,400,000 s, and 
with at least one observation no further away from the span midpoint than 35% of the whole 
span. Observations were appended to the subsegments until these criteria were met. A 
slight exception was made for data segment 17 (hourglass at ~ MJD 53,140), which only 
spans two weeks. Further observations from the end of the previous, fitted data subsegment 
were also prepended to the following subsegment when it would have otherwise spanned an 
insufficiently long time interval. An extra subsegment was generated at the end of data 
segments 7, 19, and 20 by sharing the last observation of the (otherwise) last subsegment. 

Figure 3 shows the same strong trend of v recovery (the rapid change in z>, seen immedi- 
ately after a glitch, which diminishes with time) as has been sketched in Fig. 4 of M04. Even 
in Fig. 2, the recovery of v, in the long data section following the first glitch, is apparent 
as the decrease in the negative slope of the line representing the timing solution. The long 
term trend, or steepening, from one line (timing solution frequency history) to those at later 
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times, as first reported by M04, is also apparent in Fig. 2 by glitch 15. 

Any slope, including those between the subsegment i>’s in Fig. 3, corresponds to a 
particular braking index n, as shown by the several dashed lines. The slopes are in general 
very high at the beginnings of the data segments, and in general diminish somewhat by the 
end of the data segments. The two parallel dotted lines are separated by 0.15 pHz s _1 and 
act as a guide for the changes of u within the data segments and the long term trend across 
all the data segments, which, at -0.15 pHz s -1 per five years, is an order of magnitude smaller ; 
(and opposite in sign) than the typical recovery trend of v. Although the trends within the 
segments do not stay within the bounds delimited by the dotted lines, the spans of the z>’s 
within each segment are approximately equal. This will be discussed in further detail in 
section 3.6 below. 


3.4. The Timing Solution Phase Residuals 

The phase residuals for the timing fits to all of the data segments are shown in Figure 4, 
plotted as in M04, except that we have included the phase recoveries of the third and fourth 
data segments, and most of recovery of the seventh data segment. Doing this assures that 
the fits reflect the timing behavior of the pulsar during the segment as accurately as possible, 
and are best for measuring the changes in Ay and Ai> across the glitches (see §§3.5-3. 6). As 
far we or anyone else can determine, Such post-(large)-glitch recoveries in pulsars affect the 
timing in a continuous way, unlike the pre-glitch activity such as that seen for data segments 
7, 12, 19, and 21, which will be discussed below. We have not plotted the 8 points 1 for which . 
there is no timing solution. 

Other timing irregularities are also apparent in a half dozen other data segments shown 
in Fig. 4, with both pre-glitch activity and a post-glitch recovery for data segment number 
7 (circles) - the last and incomplete segment of M04 - and in six cases have been omitted 
from the timing fits, and so are flagged by filling in the hollow characters. Segment 7 is 
plotted with greater time resolution in Figure 5 (lower). Conservatively, at least the first 
point (filled circle) is part of a recovery from glitch 6, which, at 0.46 ppm (M04), is still our 
second largest, and thus would have been expected if post- and pre-glitch activity scales with 
the size of the glitch. Data segment six contains an unphased observation as its last point, • 
supporting this hypothesis. The points beyond MJD 52,150 are shown in the inset frame of 
Fig. 5. 


r Note that one unphased point was plotted in both data segments 5 and 6, in Fig. 3 of MQ4 - the point 
belongs to data segment 5 as can be seen in Fig 2 herein. 
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In contrast, the microglitches at the end of data segments 12 and 21, plotted in Fig. 4 
as filled squares and diamonds, are very well fit by Ay’s of 0.18 and 0.29 /jHz, or gains in 
frequency of about three and five parts per billion (ppb). The duration of these features 
lasted for at least 16 and 11, and possibly as long as 30 and 18 days. This last interval 
continued, starting with an unphaseable observation at MJD 53,695 (the last unphased 
point at the end of segment 21 in Fig. 2), followed by a glitch before the beginning of the 
next segment at MJD 53,702. Both data segments 14 (Fig. 4: bowties) and 19 (elongated 
diamonds) have last points which are several sigma above the trend. Gaps of six and three 
days, respectively, follow these points. 

In addition, segment 22 is also plotted in Fig. 5 (upper), and shows evidence for quasi- 
stable behavior in the form of a steadily increasing rate of “creep” in phase, broken by a 
downward glitch, before resuming an upward trend. The timing residuals for data segment 
20 show more of a pure creep without the downward point (solid pentagons in Fig. 4). The 
creep can be considered as part of an exponentially increasing trend (see the Appendix), or 
as a slow creep of neutron superfluid vortices in the crust (Link, Epstein, &; Baym 1993), or 
both. The downward glitch behavior is rare, but a similar event apparently occurred within 
data segment 1, removed in time from any glitch (Fig. 4, leftmost squares). It has now also 
occurred in data segment 23 (Fig. 4, rightmost squares). We might speculate that the next 
glitch, expected near 2006 August 7.8 UT, is going to be a large one. 


3.5. The Glitch Amplitude-Time to Next glitch Correlation 

The glitches visible in Fig. 2 appear to have post-glitch, non-glitching time intervals 
which are highly correlated to their amplitudes. In order to test this correlation as accurately 
as possible, reliable measures of the glitch amplitudes and the time between glitches are 
desirable. We have chosen to use the frequencies of the two observation points just preceding 
and following the glitch, as evaluated/extrapolated by /from each of the two timing solutions 
from the whole data segments which preceded and followed the glitch, to determine the Ay 
across the glitch. This procedure yields two pairs of frequencies, one pair for the pre-glitch 
point from each the two timing solutions, which yields the first estimate of Ay, and then 
another similar pair for the post-glitch point, which yields the second estimate of Ay, and 
the two estimates can be compared as a check. The agreement in all 22 cases is within a 
small fraction of one /zHz. We then set the uncertainty on the time to the next glitch as a 
quarter of the sum of the gaps preceding and following the data segment between this glitch 
and the following. Thus half the entire time range possible for a glitch is given by ± one 
error, and the whole range by twice this. 
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Table 4 lists the glitch amplitudes determined as described above, while Figure 6 plots 
these amplitudes against the time interval to the next glitch, 2 for which the correlation coef- 
ficient is 0.94. The coefficient against the post-glitch stable time interval, i.e. the maximum 
time interval following the glitch when no glitches of any size are detected in the extrapo- 
lation of the timing solution (see the time interval correction in column seven of Table 4), 
is a stunning 0.96 due to improvements to the fit of glitches 2, 3, 5, 6, and 8 (these mostly 
because of the shift in glitch 1), in addition to 11, 12, 15, 17, 19, and 20, with all of these 
points plotted in Figure 7. Only glitches 4, 13, 14, 18, and 21 were worse with the correc- 
tions, but only by small amounts each. By contrast, the scatter of the amplitudes against 
the time to the previous glitch, with a correlation coefficient of -0.007, is utterly random. 

Figure 8 shows the integral of the glitch amplitude with time. The upper frame plots 
the actual - predicted glitch times. The error bars are derived from the uncertainty in glitch 
onset times and, on the left hand vertical axis, are converted into an error in //Hz through 
the slope of the oblique line drawn from the bottoms of glitch 1 to glitch 21 (but which was 
deliberately plotted offset by +10 days to allow easier comparison). The bottom. frame of the 
figure shows occasional early glitching, such as the onset of glitch 3, which may persist early 
for several more glitches, evident as the match in slope between the lower corners and the 
oblique line. Eventually an episode, during which only small (ppb) glitches occur, extends 
the time range of the data segment prior to the next large glitch, thereby resetting the glitch 
clock mechanism closer to what it was prior to the previous (early) glitch. Data segments 7, 
12, and 22, which follow glitches 6, 11, and 21, all of which end with quasi-stable behavior, 
appear to be such glitch-resetting segments. Considerable variation about this clock must 
still be possible, as glitches 13-17 alternate late and early onsets. 


3.6. The Effect of Glitching on v 

The changes in z> across the glitches, A v, and their errors, 5{Av), were determined 
in a fashion similar to those for Av, i.e., by forward and backward extrapolation of the 
timing solutions, and by evaluating v for the two observations bracketing the glitch in time. 
However, because the error in v is relatively much larger than that for u, and depends 
strongly on the epoch(s) chosen for the timing solutions, it is necessary to get Az> directly 
from the timing solutions, with epochs set to the times of the pair of observations which 
bracket each glitch (changing fit epochs was not necessary to measure the Ay’s). The errors 


2 We exclude From consideration the very small (ppb in frequency) timing irregularities at the ends of the 
data segments shown in Figs. 4 and 5 
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for Az> listed in Table 4 are taken to be minimum of the two error estimates from each of 
these pairs. 

The pattern of the glitch recovery in v shown in Fig. 3 indicates that the longer the 
interval which precedes any given glitch, the smaller the values of \v\ and thus the larger 
the possible |Az>| across the glitch. However, as mentioned in §3.3, the post-glitch \v\ does 
not seem to be able to differ from its pre-glitch counterpart by an indefinitely large amount. 
These two facts suggest that there may be a correlation between the change in i> across the 
glitch and the time interval from the previous glitch, and that it might saturate for high 
values. 

The Az>’s listed in Table 4 (and plotted in Figure 9 against the duration of the previous 
intervals for 20 of the 22 glitches) are taken as the inverse-squared error mean of the two 
values of each pair, whose difference also serves as a check of this procedure. The largest 
difference between these two is 50 ppm (0.05 on the abscissa of Fig. 9) for glitch 11. The 
next largest difference, for glitches 10 and 16, is only 0.035 in Fig. 9, and all other differences 
are even smaller. 

For glitches with 1, 000Az>z>~ 1 < 0.5 (4, 5, 8, 10, 11, 13, & 17) the correlation coefficient 
is -0.88, while that for all glitches (2-22, excluding glitch 6 due to its large error) is only 
-0.77. By contrast, the correlation of Az>z> -1 with the interval to the following glitch is 0.12 
(i.e., a negative correlation against the predominantly negative Az>’s - the same as with 
the correlation of 0.12 for Ai> vs Ax'). However, as suspected, the truly remarkable result 
is the hard upper limit of 0.075% or 0.15 pHz s" 1 for the change. The glitches near this 
limit, 2, 3, 7, 9, 15, 16, 18—2 1, and 22, all have relatively long previous intervals, bearing out 
the suspected correlation, as also does glitch 12. The correlation indicates that the interval 
which preceded the first glitch may not have been more than twice as long (~120 days) as 
the data segment for which observations have been taken, since the Az> across glitch 1 is 
0.045%. It may also indicate that the true size of Az>z> -1 for glitch 14 was likely at the lower 
end of its error range, near 0.4, rather than near the symbol plotted at 0.56 or higher, in 
Fig. 9. Noticeable improvements (a few to several units or so) in the chi-squares of the runs 
achieved, e.g., by using the polynomial fit parameters to initialize a minimum sum-squares 
routine, will typically result in changes to A v in Table 4 by only 1-2x10" 14 Hz s -1 . However, 
in none of these cases did the Az> ever exceed 0.15 pHz s" 1 . 



4. Recapitulation 


The observations of J0537 have revealed several unique and remarkable features “in 
need of explanation.” The first is the tight correlation between the size of the glitches (A y) 
with the time interval to the following glitch, with a slope of about 400 days ppm -1 . The 
correlation with the post-glitch stable time interval is even tighter. Second, the microglitches 
which precede the large glitches can be due to a simple small increment to the pulse fre- 
quency, y, or an initially gradual, but more than linear change in frequency with time, or 
a combination of both. The data are also consistent with each glitch being followed by a 
smooth, but relatively small recovery involving the decay of only a few percent of the initial 
Ay. 

Third, the long term trend of \v\ with time is increasing , and thus the timing age of 
J0537, -0.5yy -1 , continues to decrease at a rate of nearly one year per year. Another way 
of stating this is that the gain in |yj across the glitch, typically about 0.15 pHz s -1 , unlike 
the behavior of z> in Vela noted by Alpar (1998), is not completely given back before the 
following glitch, with 10 percent of the gain remaining. Fourth, the change in y, or Az>, is 
correlated to the time interval to the preceding glitch, but clearly saturates at 0.15 pHz s -1 . 
And finally, the RXTE-band (~2-20 keV) pulse profile has only a single pulse with a flat 
maximum for at least 0.02 cycles, or 320 ps. 


5. Discussion 

In §5.1 we discuss the strong correlation between glitch size and the time to the next 
glitch, and explain how this helps in discriminating between current models applied to other 
glitching pulsars. To this we add a discussion of the change in spindown rate, Ay, across 
and between the glitches in §5.2.1, and a brief discussion on glitch latency intervals in young 
pulsars in §5.2.2. The remaining longterm increase of |y|, i.e., the increasing spindown with 
age, is then discussed in §5.2.3, including the possibility that the magnetic pole of J0537 is 
migrating away from its rotation axis, and the necessity for a glitch latency interval in the 
youth of J0537 similar to that observed in the present day young pulsars. Finally we finish 
this section with a brief discussion of the pulse profile, in §5.3. 


5.1. The Glitch Size-Time to next Glitch Correlation 

It is generally accepted by now that the moment of inertia contained by the neutron su- 
perfluid within the crust of an NS amounts to only about 1% of the total, which also includes 
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both that of the superfluid interior and the solid crust. Some glitches, including those of 
the Crab pulsar and B0540-69, may be “crustquakes,” where the equilibrium configuration 
(EC) for the solid crust departs from its geometrical configuration as the pulsar spins down 
until eventually the crust cracks and settles. 

Of course, settling of the crust alone can not possibly be responsible for all glitches in 
Vela, as the NS would run out of its supply of rotational oblateness after a few hundred 
years. Also, sufficient crust settling to produce a frequency change of only a fraction of 
a ppm, for a glitch in a low luminosity pulsar such as Vela, would involve the release of 
an easily observable amount of energy (likely in the soft X-ray band), which was not seen 
(Andersen & Ogelman 1995; Seward et al. 2000), In addition to the glitch models mentioned 
in this work, several alternative models for glitching in pulsars exist (see, e.g., the discussion 
and references in §6 of Dall’Osso et al. 2003). Other models have been proposed more 
recently (Charnel Sc Carter 200.6; Jahan-Miri 2006a, b; Negi 2006; Peng, Luo, Sc Chou 2006; 
Zhou et al. 2004). 

In glitching pulsars, the solid crust of the NS spins down continuously in between 
glitches, as the magnetic field is locked to it. 3 * During this time at least part of the neu- 
tron superfluid within the crust is decoupled from its spindown, and thus becomes a second 
(and separate) stream. Thus, as the solid crust slows down, the vortex density of the sec- 
ond superfluid stream exceeds what would be appropriate for the solid crust (and the first 
superfluid stream) by an ever larger amount with time. 

Andersson et al. (2004) have suggested that such a two-stream superfluid system is 
subject to an instability, analogous to the Kelvin-Helmholz instability, that acts when the 
difference in effective rotation rates between the two streams is greater than a certain amount, 
Oja - . During a glitch this causes a transfer of at least part of the extra angular momentum in 
the second, and faster rotating superfluid stream, to the solid crust and the first superfluid 
stream within it, spinning up the pulsar by a certain rotational difference, A v. Another 
view holds that an excess of vortices from “traps” located within the crust, transfer angular 
momentum to it when triggered by a variety of potential mechanisms, including crustquakes 
in addition to an (see, e.g., Alpar. 1998; Alpar et al. 1993; Epstein Sc Link 2000). We 
explore this issue further below. 

The actual mean size of a glitch in J0537 after about 10 7 s (115.7 days) would only be 
about 18 pHz, or just over one day of spindown (see glitch 3 of Table 4), which is nearly 
exactly one hundred times smaller than the accumulated spindown over this time interval, 


3 Even now, it is not yet clear how or if the interior superfluid of the NS is coupled to the solid crust 

(Sedrakian & Sedrakian 1995; Andersson, Comer, & Prix 2004). 
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and thus supports the hypothesis of angular momentum transfer from an NS component 
with 1% of the moment of inertia (using 399.1 days per glitching ppm for J0537 yields 0.9% 
of the spindown is being reversed by the glitching). For the Vela pulsar, 4 the glitches average 
2.0 ppm, or 22 /uHz, very close to the average glitch size for J0537, but this amounts to 17 
days of spindown, or just over 1.5% of the accumulated spindown over the 2.8-year mean 
interval between its glitches, which again- supports angular momentum transfer from a 1% 
NS component. The difference between 1. 5-1.7% in Vela and 1% in J0537 may be due to 
the difference in their Q values, 5 as Wang et al. (2000) report 0.38 for the 1996 October 13th 
glitch. 

The coincidence between the mean absolute glitch sizes for J0537 and Vela suggests 
that the neutron superfluid in the crusts of all NSs, 6 when triggered by some mechanism, 
which may not be the same for all pulsars, favor dumping some of, or even more than (in 
the case of finite Q), Ua-. Relatively “large” glitches, amounting to gains of 32.4 and 28.0 
fj. Hz, still essentially within the range observed for both Vela and J0537, have also been 
observed from the 2 Hz PSR J1806-2125 (Hobbs et al. 2002) and the 4 Hz PSR J1614-5048 
(Wang et al. 2000), respectively. Although these increases in spin frequency are the same 
to within a factor of 2, the corresponding increases in energy span more than an order of 
magnitude, which lends yet more support for the storage and exchange of angular momentum 
by vortices in the neutron superfluid which is thought to lie within the crust. 

Although glitch 1 of J0537, at 42.2 /tHz, is the largest glitch ever observed, the 34.5 and 
34.3 pHz Vela glitches of 1978, July 13, and 2000, Jan. 16 come in a close second and third, 
followed by the 32.4 /jHz glitch from PSR J1806-2125 as a very close fourth, and still ahead 
of all of the other glitches seen from J0537 to date. All known glitches with Av > 5 pHz, 
are listed 7 in Table 5, and histogrammed in Figure 10, which shows a clear maximum size 
for all glitches of 40 pHz. This lends further support for an u v, between the solid crust 
and part of the superfluid contained within it, and this maximum can not be much bigger 
than 40/uHzxl00, or 4 mHz, i.e., one “revolution” every 250 s (Anderson & Itoh 1975). This 
would correspond to 8 years of spindown for Vela, but only about 8 months for J0537. 


4 Lyne et al. [2000] report 1.7% for most pulsars. 

5 Q is the fraction of the Au gained across a glitch which decays afterward. 

6 The existence of a maximum in Au for the crustal neutron superfluid also implies that the glitch of AXS 
J161730-505505 reported by Torii et al. (2000) likely occurred closer to MJD 50,600 and had an amplitude 
near 26.2 jullz, instead of between MJD 49,300 and MJD 50,000, when its amplitude would have had to have 
been between 42 and 115 /iHz. We note that the timing age of AXS J1 61 730-505505 is 16,000 years. 

7 Zavlin, Pavlov, and Sanwal (2004) report three possible glitches in the 424 ms pulsar, 1E1207-5209, all 
with likely At/s > 5 /tHz, but other interpretations for the timing irregularities also exist. 
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The NS crust in Vela is nowhere near as stressed even after its 2.8-year mean inter- 
(large)gliteh time interval than it is in J0537 after only four months. Specifically, the crust 
of Vela would have had to mostly settle by only 0.007 cm after 2.8 years, whereas the crust 
of J0537 would have had to mostly settle 0.06 cm in four months, as the settling distance, 
JR, where R is the radius of the neutron star and is assumed to be 12 km, 8 scales as vi>. The 
factor of 70 difference in crust settling rates between the two may mean that less than one 
large Vela glitch in 15 is initiated by crust cracking. After J0537 and the Crab, the glitching 
pulsar with the next highest i>v (down by a factor of 18) is the 15 Hz J0205+6449. The 
pulsar with the lowest vv (down by a factor of 12,000) is the 2 Hz J1806-2125, but its v is 
down only by a factor of 200. Clearly most glitches in most pulsars can not be triggered by 
crust cracking; many of those in the Crab and B0540, and almost all of those in J0537 likely 
are. 9 

Since the large glitches in J0537 are probably initiated by crust cracking events, the 
extreme linearity in its glitch size-post glitch quiet timing interval could be 4ue to the 
steady divergence of rotation rates, and of the two superfluid streams, and crust 
cracking occurring much more frequently than its large glitches, 10 which only triggers them 
when v 2 — v x > w CT . The alternative, of having the amount of angular momentum dumped 
from many vortex traps proportional to the amount of settling in its solid crust, is harder to 
justify, as many vortex traps with differing saturations would result in a chaotic glitch- time 
relationship. The crust settling in J0537 would have to be a global, rather than a local event, 
dumping the vast majority of traps to the same level, which seems unlikely. 

Another alternative would be to have just one trap (Franco et al. 2000; Jones 2002), but 
the pre- (large) glitch microglitches visible in Figs. 4 and 5 contraindicate this alternative. The 


8 The mean settling rate [drop/time] goes as 4 tt 2 z>vR 4 (GM) -1 , and thus would double for R =14 km. 

9 J0537 is so bright in most of the X-ray band, however, that the excess luminosity due to its crust settling 
might easily escape notice. Assuming, in the case of J0537, some 10 7 s worth of spindown causes the settling 
of a 1.2x10® cm radius by 0.054 cm and a glitch of 17.8 /zHz, the specific energy gained from the gravitational 
field of 1.3xl0 14 cm s -2 would be 7xl0 12 ergs/gm. Using a specific heat estimate equal to that of water, 
or 4.17xl0 7 ergs gm -1 K -1 , gives a temperature rise of 14.5 eV. The temperature rise for J0537’s biggest 
glitch would be 34.5 eV. With an NS surface temperature of 80 eV, these rises might not be detectable, 
unless the specific heat is overestimated by a large factor. Link and Epstein (1996) have argued that such a 
medium deposition of energy within the NS crust could have caused the large glitch seen in Vela during 1988, 
Dec. 24 (Lyne, Smith, & Pritchard 1992), and have also suggested the same mechanism for Crab glitches, 
the difference in magnitude being due to the temperatures of the two NS crusts, i.e., age. One large Vela 
glitch out of the dozen known is consistent with the >15:1 ratio estimated above. Unfortunately, no X-ray 
observations bracket the 1988 Vela glitch. 

10 The cracking event glitch rate in Vela then would be larger than 1 in 15. 
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relatively poor linearity observed in Vela glitches would then result from two-stream rota- 
tional differences significantly greater than w cr before cracking could trigger a glitch, leakage 
due to vortex currents in the inter-glitch interval, or both. We note that a more complicated 
picture of the glitch processes within Vela is asymptotic to this linear relationship as the 
Ai>’s across its glitches become equal, which they very nearly are for the long term recov- 
ery components of the larger glitches of Vela (Alpar et al. 1993), certainly more so than for 
J0537. 

Although it has been suggested that cracks play no role whatsoever in glitches (Jones 
2003), the arguments are based on the difficulties of producing open cracks within the crust, . 
while closed cracks are perfectly free to slide and grow, given sufficient stresses. The crust 
does have a certain mechanical strength, even with cracks, but this can not literally hold 
it up as it strains against gravity for any indefinite period of time as the stresses continue 
to increase. At stresses greater than the yield stress, 11 the solid crust may have a range of 
higher stresses wherein it undergoes plastic deformation (Ruderman 1991a, b,c; Ruderman 
et al. 1998). In this stage the crust will slowly dissipate most of the gravitational energy 
involved in its settling, and thus avoid a sudden, large unobserved release of gravitational 
energy during a glitch (from whatever settling is left over). The older a pulsar is, the more 
it has glitched in the past, and thus the more extensive the network of cracks is likely to be. 

Eventually, for J0537, the stresses in the crust will exceed the stability criterion for 
a population of cracks which- is always present, and these cracks then grow quickly in an 
unstable fashion (see, e.g., Dienes 1985 and references therein) until they relieve the stresses 
and subsequently stop growing, possibly coalescing into a fault which crosses the equator 
(Franco et al. 2000), and if — iq > Wcr, triggering a large glitch. It is also possible that, 
during this process, the cracks will network throughout the entire crust as they grow and 
intersect each other. The mechanics of cracks in solid materials, including the NS solid crust, 
is discussed further in the Appendix. . 

A glitch can occur slightly early if the previous glitch left cracks which were larger than 
those which usually occur, and these become unstable earlier (at lower stress levels), but 
would not cause a glitch which was too early because — iq < In these cases the stress 
levels at which crack growth become unstable will be lower. A glitch could also occur later 
if the population of larger cracks is relatively depleted. 


11 For stresses below the yield stress, a solid body undergoes a strain which is strictly proportional to 
the applied stress, and this strain vanishes when the stress is removed. This is Hooke’s Law for the elastic 
behavior of solids. At higher stresses the extra strain in the solid is permanent, and this is due to plasticity, 
crack growth, or a combination of both. Another strain component which can vanish when stresses are 
removed is due to crack opening. 
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5.2. The v Behavior 
5.2.1. The A v Behavior 

We are still faced with the A v across the glitches, its correlation with the interval prior 
to the glitch, the apparent saturation value of 0.15 pHz s -1 , and its long term increase in 
absolute value, in addition to the behavior of v between the glitches. If we consider that 
J0537, at 4,000 years of age, is a middle-aged pulsar which is much older than the Crab, and 
yet still less than half as old as Vela, the observed behavior of at least the v is what might 
be expected, at least empirically. 

The glitches in the Crab pulsar (and in B0540) are sufficiently small and infrequent,, 
with a dozen glitches with Aw' 1 from 4xl0~ 9 to 6xl0~ 8 (Wong et al. 2000) between 1969 
and 1999, 12 that they could only be due to a few vortex traps within the crust. Thus the 
excess angular momentum dumped into the crust when a glitch occurs in the Crab is small, 
and the glitch size is small (Alpar 1998; Alpar et al. 1994,1996). 

In addition much of the gain in Avis' 1 across its glitches, up to 4xl0 -4 , persists indefi- 
nitely, initially prompting some (Allen & Horvath 1997; Link & Epstein 1997; Link, Franco, 
& Epstein 1998; Ruderman et al. 1998) to attribute this to a gain in the magnetic dipole 
moment due to an increase in angle between the magnetic and rotation axes (see §5.2.3). 
However, the pulse profile of the Crab already has two peaks, and is thus hard to recon- 
cile with pole migration. Perhaps a more likely possibility, as first suggested by Alpar et 
al. (1996), is that the initial crackup permanently destroys the ability of part of the Crab’s 
solid crust to store vortices (creating “vortex depletion regions”) thus reducing its effective 
moment of inertia, which, for a constant braking torque, leads to a permanent increase in 
spindown. On the other hand, Ruderman (2005) attributes this effect to a permanent change 
of braking torque across the glitch due to an increase of the component of the dipole moment 
perpendicular to the rotation axis, which, in turn, is caused by migration of magnetic flux 
tubes and vortices, rather than to any physical migration of the crust along the NS surface.. 

In the much older (wrt the Crab) Vela pulsar, the NS solid crust may be saturated 


12 Although the 1969 glitch of the Crab (MJD 40,493 - Boynton et al. 1969) is generally quoted as 4xl0 -9 
in Aw' 1 , the observing coverage very near its time of occurrence was sparse, and there is good reason to 
believe that the actual Aw' 1 was a full 1.0 . This is because a glitch was observed, in data taken by 
Rem Stone for J. E. Nelson and J. Middleditch, around 1 August, 1971 (very early in the [optical] observing 
season, which started earlier yet on 24 July), which was very accurately measured to be 4.2xl0~ 9 in Avv~ l . 
This glitch overlaid the 1969 glitch nearly perfectly when plotted on the same time scale but a different phase 
residual scale by a ratio of exactly 2:5. The timescale of an exponential decay of Av for this glitch, r c , was 
2.4 days, in good agreement with r c = 3 days reported by Wong et al. (2001) for its weaker glitches. 
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with cracks, allowing a considerable number of traps, and/or a substantial second superfluid 
stream, so that the excess angular momentum dumped into the solid crust during a glitch 
spins the whole star up by an average near 22 fj , Hz every 2.8 years. The extra Aui>~ 1 gained 
across its glitches, 4-7xl0~ 3 , appears to completely die away in a what seems to be a linear 
fashion (i.e. with constant braking index) through vortex currents (Alpar 1998), and the 
subsequent glitch is primed to occur when this process finishes. For Vela this means there’s 
also a strong correlation between Az> across a glitch and the time to the next glitch, in 
addition to the same correlation for Av (Alpar et al. 1993). 

The situation for J0537 differs either because the crust frequently doesn’t wait for cur- 
rents to finish giving away all of the gain in j Az/j before causing a glitch, or because of an 
increase in the strength of the magnetic field, or a small, permanent gain caused by realign- 
ment of the NS magnetic moment, or any combination of these. This does not require that 
the vortex currents in J0537 and Vela differ in any fundamental way, as the jAz/j gained 
from the glitch in J0537 subsequent to the early “giveaway” glitch is usually correspondingly 
more modest, as can be seen in Fig. 3. 

Although the glitches of the adolescent pulsar, J0537, average about six times smaller 
than those of Vela in Aw" 1 , they are an order of magnitude more frequent, so that its 
glitch activity parameter, 13 is 0.9 ppm yr _1 as compared to Vela’s 0.7 ppm yr _1 , and at 
0.1-0.7xl0~ 6 , they are still larger than those of the Crab by about one order of magnitude. 
About 90% of its gain in Az>z> -1 (as opposed to 100% in Vela), of up to 7.5xl0 -4 (about an 
order of magnitude smaller than the relative gains in Vela) “decays away” prior to the next 
glitch, the only other difference from Vela being the softening of the J0537 braking index 
between its glitches. Since the absolute size of the glitches in J0537 are just as large, and in 
one case larger than those in Vela, the crust of J0537 is likely just as saturated with cracks 
as Vela. This may be more of a reflection of enhanced crust settling in J0537 due to its 
history of higher spin rates, than it is of actual age. 

However, the relative amount of the gain in Avv~ l across a glitch in Vela that decays 
away afterward, or Q, at 0.38 (Wang et al. 2000) is about an order of magnitude larger than 
in J0537. The integrated extra spindown due to a typical glitch of J0537 of size £ ppm, with 
a linear Av recovery of 0.21 £ pHz s -1 over the 400 £ days until the next glitch, would slowly 
remove 0.059 £ 2 ppm from the gain in v across the glitch, during the interval to the next 
glitch, which is consistent with the barely visible recovery following glitch 1 shown in Fig. 2. 

This decay may be caused by vortex currents from the interior superfluid to that in 


13 The glitch activity parameter is defined as: A g = Y,Av g /v/At, where £A^ 3 is the sum of all spin 
frequency gains to v from glitches, over a timespan, At. 
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the crust, during the interval following a glitch, which makes up for an excess of vertices 
dumped from traps and/or a second superfluid during the glitch. Interestingly, the agreement 
between the maximum values of the A i> of 0.15 pHz s _1 involved in the long term recovery 
following the glitches of Vela, J0537, and the Crab (the 1989 July 13 glitch - Wong, Backer, 
& Lyne 2001) is nearly perfect, even better than that of the maximum values of Av among 
all pulsars. This is what would be expected from one NS to the next, as crusts, vortex 
currents and moments of inertia are basically the same, perturbed possibly only by spin 
rate, magnetic field, and differences in mass and temperature. 


5.2.2. Young Pulsars and the Latent Interval 

Livingstone, Kaspi, and Gavriil (2005) have called attention to the two orders magnitude 
difference in the glitch activity parameter, between the Crab pulsar and B0540.as evidence for 
significant differences in the internal structure of neutron stars. Although the most obvious 
difference between these two is that between their pulse profiles, the harmonic content of 
B0540 (Seward, Hamden; &c Helfand 1984; Middleditch & Pennypacker 1985) is similar to 
the optical pulse profile of Vela (Wallace et al. 1977), and distinct inclinations for the two 
sources could easily account for the rest. It is more likely that young pulsars such as the Crab, 
B0540, and B1509-58 don’t glitch very frequently at first, having to rely on settling of their 
crusts in order to form their initial cracks (generating vortex depletion regions and vortex 
traps and/or more second superfluid domain in the process), which then help in the further 
crackup of the NS crust. In this case, the crustquake rate would increase exponentially (at 
least initially), with the product, vv, and the age of the pulsar, r: 

vAg ~ exp (t/T vv/{vv)crab),_ (4) 

where T is an e-folding timescale of glitch activity. If we use the timing age of the pulsar, 
-0.5vz> -1 , in place of r, then equation 4 becomes: 

vA g oc exp (-0y/(w)c ril i/T), (5) 

and the e-folding time for B0540 would be three times as long as that of the Crab, due to its 
product, vv, being three times smaller than that of the Crab, plus 25% because the timing 
age, r, of B0540 (1,550 years) is 25% longer than that of the Crab (1,240 years). 

The absolute glitch activity, EAv g f At, for B0540 is really 75 times less 14 than that of 
the Crab, but this can be made up if T = 170 years. For B1509, with a 6.6 Hz spin, a 


14 The Crab’s vA g is taken from Livingstone et al. (2005) to be 0.3 n Hz yr x , whereas that of B0540 is 
0.004 pHz yr -1 . 
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6.7xl0 -11 Hz s -1 spindown, and a 1,670 year timing age, its vA g is only 28% of that of 
B0540, consistent with no glitches having been observed to date from B1509, both pulsars 
having been discovered in the early 1980s, If no starquake is observed from B1509 in the 
next few decades, then perhaps equation 4 will need a pre-exponential factor of z>, which is 
certainly needed for the glitching rate of mature pulsars like J0537 and Vela. 

This is a simpler model than that employed by Alpar and Baykal (1994). The starquake 
rate must turn over and flatten out if the vA g predicted for J0537 is to match the actual 
value of 56 //Hz yr _1 , i.e., only a factor of 187 higher than that of the Crab. This vA g is 
only a factor of 7 more than that of Vela, as compared to the factor of 12.5 expected from 
the ratio of their z>’s, both pulsars being old enough so that their cracking growth functions 
may have reached, the same high constant value. This discrepancy is not serious and can be 
completely resolved if the Q in Vela is about 0.46, a value close to that measured, by some, 
or if Vela is old enough to have a higher value of the cracking growth function than that of 
J0537, or a combination of both. •• ; • 


5.2.3. Magnetic Pole Migration 

Three facts support J0537 as being born, and possibly still remaining, close to an aligned 
rotator: the singly-peaked pulse profile; the low value of its v with respect to that of the 
Crab pulsar, 15 which has almost twice the i> and less than half the u; and the long term 
increase of its |z>|. The increase in the magnitude of the (negative) z> amounts to about 0.15 
pHz s _1 every five years, or 0.95 xlCT 21 Hz s~ 2 , but is still uncertain by about 25%, Settling 
of the NS over this time period due to spindown will produce a 0.1% reduction in the 8.4 
m NS equatorial, centrifugal bulge, or 0.84 cm, and a change in the dipole moment due to 
scale size changes a thousand times smaller, if any at all because the magnetic axis is not 
close to either the rotational equator or pole. 

If the entire stellar moment of inertia can drop by 0.1% over five years, then, provided 
the braking torque doesn’t change, this could account for the 0.075% increase in \v\ in the 
same period. However, if the effective magnetic dipole moment of the NS also drops, the 


15 Although the two pulsars have identical spindown luminosities (pltPlvv) of 5xl0 38 ergs s -1 , close to 
the Eddington luminosity for an NS, SN 1986J, which is now producing 200 times the luminosity of the 
Crab Nebula in the 2 cm band (Bietenholz, Bartel, & Rupen 2004), likely exceeds this. in total luminosity 
by a large factor. Thus the spindown of J0537 is not smaller than that of the Crab because its spindown 
luminosity, and that of all pulsars, L s d, is not likely to be limited to L e dd, even by processes which we do 
not yet fully understand. We also note that Shukre, Manchester, and Allen (1983) have suggested that the 
SS 433 binary system contains an aligned, rotating, accreting NS. 
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reduced braking torque could cancel the effect, but if it increases, it could help produce it 
(see, e.g., Ruderman 2005). If this effect is real in J0537, then it is because its has the highest 
product, \vv\ of 1.2xl0 -8 Hz 2 s -1 . However, the same product for the Crab pulsar, l.lxlCT 8 
Hz 2 s -1 , is a close 2nd, and the Crab clearly has a robust, positive v of 1.2 xlO -20 Hz s~ 2 . 
Still, some have suggested that the intrinsic magnetic fields of pulsars may increase with time 
(see, e.g., Blandford, Applegate, & Hernquist 1983), though until now, there has been little 
evidence to support this assertion (Wakatsuki et al. 1992), and the magnetic field of the NS 
remnant of SN 1986J is clearly very strong after only 20 years (Bietenholz et al. 2004). The 
effect in J0537 is consistent with a growth in its magnetic field of 0.375 xl0~ 4 yr -1 , or 35 
MG yr -1 . 

.There remains only the migration of the magnetic axis as the possible cause of the 
long term increase of |z>). As mentioned in the previous section, some have initially argued 
that the persistent spindown gains in the Crab pulsar following its glitches are not caused 
by pole migration (see, e.g., Link, Epstein, & Baym 1992), but instead by the formation 
of new vortex depletion regions, (Alpar et al. 1996), However, this only works for young 
pulsars which haven’t glitched at the high rates that Vela and certainly J0537 have for many 
thousands of years. 

Following arguments similar to those made previously by others (Allen & Horvath 1997; 
Link & Epstein 1997), and assuming that the spindown of J0537 is proportional to the ef- 
fective magnetic moment, |m e //| = |m x v\\ v\~ l , squared, or |m e //| 2 , and that this quantity 
is in turn proportional to sin 2 a, where a is the polar angle of m e // from the rotation axis, 
and is also a linear function of time, we have: 


a(t) = at + a 0 , 

(6) 

i/ — $ sin 2 a{t) 

Vpsr 

(7) 

/? sin 2 a(r psT ) = — 2 x 1CT 10 Hz s -1 , 

(8) 


where n is the braking index (for now assumed to be 2.5 - the same as is characteristic 
for the Crab pulsar), a is the time derivative of a(t), r pST is the age of J0537, ao is the 
initial value of a(t) for time 0, or —t p „ from the present, and /3 is a negative constant with 
dimensions of v. As the magnetic pole migrates away from the rotation axis, the increase in 
the effective magnetic moment produces the increase of the spindown as a function of time 
that is evident in Fig. 3. Thus the time derivative of equation 7 is equal to the long term 
trend seen in Fig. 3: 

v = 2a/3 sin a(r pST ) cos a(r psr ) + ni/ psT /3 sin 2 a(r p3T ) /u psr = —0.95 x 10 -21 Hz s" 2 . (9) 
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We can remove the second term from the middle of equation 9 by subtracting an 
assumed known present day power law contribution, to v for J0537, with n — 2.5, of 
ni>p Sr fi sin 2 a(r p3r ) / v p3T = nv p3r /v p3T = 1.6 x 10~ 21 Hz s -2 , to get: 

i> a = 2a/3 sin a(r p3r ) cos a(r psr ) = —2.55 x 10 -21 Hz s~ 2 . (10) 

Dividing equation 8 by equation 10 to eliminate /? gives; 

tana(r p3r ) _ ? g x 1Q io a> Qr 2) 475 yearS) (11) 

a transcendental equation which can be solved for a, with a(£) as the simple linear function 
of time given in equation 6, and if ao and r p3T are known. 16 Values for a and ao which make 
a(r pjr ) exceed 90 degrees are naturally excluded. We can integrate equation 7 to extrapolate 
backward and forward to solve for v{t) and u(t). 

Since the absolute value of the tangent function is always greater than that of its argu- 
ment, and d, and ao are, barring migration toward the rotation axis, positive, equation 11 
indicates that the age of J0537, r p3T , or at least, r a < r p3T , the duration for which there is a 
constant a, can not exceed 5,000 years. For reasonable values for ao and a with r a — 4,000 
years, equation 11 is a problem to solve. For example, with ao = 0.2 radians, no : value of d 
exists that gets the left hand side of equation 11 below 3,368 years, as long as we assume r a 
= 4,000 years. 

The magnetic pole migration rate, necessary to render equation 11 solvable, is plotted 
against its initial obliquity for various values of latency times, r psr - r a for t P3T = 4,000 years, 
in Figure 11 (lower). Four vertical lines for ao = 0.1 to 0.4 radians, in steps of 0.1 radians, 
intersect the lower portions of the curves for 22 (believable) values of d. From equation 6 
we have, for the present day magnetic obliquity, a psT : 

a psr — dr a ~f" ao, (12) 

which we can rewrite for the duration of the active migrating interval, r a as: 

Ta — ao/d - b a psr / d. (16) 

Thus, given a psr and d, the relation between ao and r a is a straight line with a slope of 
— 1/a and an intercept of a p3T /a. This relation is plotted in the top frame of Fig. 11 for 


16 For a braking index, n, of only 2.0, near the lowest known established, glitch-free braking index for any 
young pulsar (Zhang et al. 2001; Cusumano, Massaro, & Mineo 2003; Livingstone et al. 2005), the right 
hand side of equation 11 increases to 2,830 years. 
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cxpsr from 10 to 60 degrees in steps of 10 degrees, and illustrates those values of a psr and cto 
which are necessary to keep the non-glitching interval believably small. 

The time histories of j'(t) and z>(t), corresponding to the pairs of (cxo,a) from the lower 
frame of Fig. 11, are plotted in Figure 12. The figure shows a small range for the spin rate, 
v, of J0537 at birth, between 75 and 80 Hz, and a much wider relative range for v. 

5.3. The Pulse Profile 

The peak of the pulse profile of J0537 appears to be flat for at least 300 fjs (Fig. 1). 
Following the line of argument from Golden et al. (2002) to test whether J0537 should be 
able to produce high energy pulsations, given a 0.02 cycle extent of the plateau, we determine 
that the magnetic field of the emitting region, B em in Gauss is given by: 

Bern — 10 10 E, (.14) 

where E is the energy of the peak of the pulsed synchrotron emission in keV., However, 
the Lorentz 7 is only 8 from the same calculation, so it is hard to know if this analysis is 
appropriate for J0537. Assuming that it is, and extending a dipole field of 0.925 tG, and 
a radius of 12 km, gives 54 km for the location of pulsed keV emission. The field strength 
at a light cylinder with a 770 km radius would then be 3.5 MG, and thus J0537 would be 
expected to produce optical pulsations. However, if J0537 is a nearly-aligned rotator, the 
actual surface magnetic field may be near 3.5 tG, or higher. In that case the pulsed keV 
emitting radius moves out to 86 km, and the magnetic field at the light cylinder increases 
to 14 MG. 

If the field does actually fall through 2-3 xlO 7 G before reaching the light cylinder, by 
the logic of some (Cheng, Ho, & Ruderman 1986a, b) , one might expect an outer gap region 
where pulsed optical emission could arise. Gil, Khechinashvili, and Melikidze (2001) model 
the pulsed optical emission from Geminga (Golden et al. 2002) assuming that, like J0537, it 
is a nearly aligned rotator. However, O’Connor, Golden, and Shearer (2005) argue that, in 
the case of J0537, its pulsed radiation is synchrotron self-absorbed at optical wavelengths, 
with my ~ 24. Unfortunately, the limits on optical pulsation are not very stringent, with 
the magnitude 23.4 limit from Mignani et al. (2000) derived from imaging alone and close 
to one solar luminosity, as compared to the Crab pulsar with 4 solar luminosities. 



-24- 


6. Conclusion 

The 62 Hz pulsar in the LMC, J0537, is unique among all others. It is the fastest spinning 
young pulsar, and the most actively glitching pulsar known, with a gain in frequency of 56 
/iHz, as compared to Vela’s 8 /rHz, yr -1 . The extreme linearity of the glitch-size/post-glitch 
stable time, in addition to other quasi-stable behavior, is consistent with the crack growth 
mechanism in the NS solid crust for the initiation of its glitches, and a two superfluid stream 
instability within the solid crust, as the cause of the variation in glitch magnitude. Crack 
growth mechanisms may also cause glitching activity in young pulsars, such as the Crab and 
B0540, to increase exponentially with a 170-year timescale and the product, vv[vvc Ta b)~ l , 
and this may be detectable in the Crab after a few more decades. Among all pulsars, the 
important .quantity involved in glitches is angular momentum transfer, i.e., Av, rather than 
Ai/j/ - 1 or vAv, and 40 uHz appears to be an absolute maximum for this quantity. This 
represents an w^., or a difference in rotation rate of one cycle per 250 s between the NS 
solid crust, and the superfluid vortices within it. The moment of inertia of the crustal 
superfluid is close 1% for J0537, lower than that ascribed to Vela and a few other pulsars, 
and the difference may explain why Q, the fraction of Av gained in the glitch which decays 
afterward, is an order of magnitude higher in these others. 

The 62 Hz rotation frequency of J0537 is over twice that of the 30 Hz rate of the Crab 
pulsar, but yet its spindown rate is only slightly greater that half that of the Crab, and just 
barely larger than that of the 20 Hz rotator, PSR B0540-69. Clearly its effective magnetic 
dipole moment is only about 25% those of Vela, the Crab, and B0540. It may just have a 
weaker magnetic field than the others, but its singly-peaked pulse profile argues that, when 
it was born 4,000 years ago, it was as a nearly-aligned rotator, spinning at a rate between 
75 and 80 Hz. 

The longterm increase in the magnitude of its (negative) spindown, \v\, near Iff" 21 Hz 
s~ 2 (or 2.55 xlO -21 Hz s -2 after correcting for braking by magnetic dipole radiation and/or 
pulsar wind processes), supports the interpretation of magnetic pole migration away from 
the rotation axis over time, by about one radian every 10,000 years, or one meter on the NS 
surface per year, or even less. This rate agrees with the estimate, made in the Appendix, of 
the crack growth rate necessary to trigger the glitches in J0537. Like the Crab, or B0540, 
which, after 950 and 1,550 years, are both still in a low glitch-activity state, the interval of 
relatively low glitch activity in J0537 may have lasted 1,500 years, although realistic values 
for the initial obliquity and migration rate also exist for an interval as short as 400 years. 
The spindown rate, \i/\, for J0537 at birth was likely at least 30% less than it is today. 

There also appears to be a maximum absolute value for the (almost always negative) 
change in spindown, |Az>|, across (and for J0537 and Vela also equal to that involved in the 
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longterm recovery following), a glitch, in all pulsars, of 0.15 pHz s -1 (i.e., a minimum value 
of -0.15 pHz s -1 ). This is equivalent to a maximum value for vortex currents within pulsars, 
though exactly how or why this occurs is still unknown. Finally, the single peak of the J0537 
pulse profile is flat at its top for at least 0.02 cycles, or 320 /j,s. 

J.M. thanks Drs. Ali Alpar, John Dienes, Aaron Golden, Naoki Itoh, David Pines, and 
Mai Ruderman for helpful, guiding conversations, and Tom Bowles of Science, Technology, 
and Engineering, and the Institute for Nuclear and Particle Astrophysics and Cosmology, for 
funding. E.V.G.’s research is supported by NASA grant NAG05GR11G. This research has 
made use of data obtained through the High Energy Astrophysics Science Archive Research 
Center Online Service, provided by the NASA/Goddard Space Flight Center. This research 
was performed in part under the auspices of the Department of Energy. 
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7. Appendix 

From Dienes (1983) as corrected by Rice (1984), a penny-shaped crack will become 
unstable to rapid growth if the following inequality holds: 17 

(cr -r) 2 c > (1) 

where e is the crack radius, a and r are the traction and interfacial friction across the crack, 
v p is Poisson’s ratio 18 , j $e is the specific surface energy (in, e.g., ergs cm -2 , or dynes cm -1 ), 19 
and fi is the shear modulus. 20 As the crack radius, c, increases, unstable growth will occur 
for lower stress levels. 

Melting and “healing” of the cracks due to interfacial friction may also occur as their two 
surfaces slide against each other. Others may stop growing when they intersect neighboring 
cracks at a high angle (“T” cracks). Still others may rotate or translate within the solid 
crust so that the traction stress is reduced across the largest and hence, most unstable cracks, 
and others yet may disappear through subduction below the crust (Ruderman 1991a). As 
occurs with faults in the Earth’s crust (see, e.g., Scholz 2002), the crack intersections in the 
NS crust and other local variations prevent the crustquake process from relieving all strains 
necessary for it to settle to its current EC. 

As stated in the main text, data segments 7, 12, and 21, which may restore the “original” 
glitch timing shown in Fig. 8, all have a period of quasi-stable behavior during this interval, 
instead of a prompt (large) glitch. This behavior can be physically interpreted as the NS still 
having the “guilty” cracks which caused the onset of the previous glitch, except that these 
have been partially “repaired” by melting or some other mode of fusion for closed cracks. 
The repair process leaves the crack radius smaller, thus revising the failure stress Upward 
toward the nominal value. As mentioned above, another possibility is rotation or translation 
of the crack over several consecutive glitch intervals until the traction stress at the time of 
early onset is no longer sufficient to cause the crack to grow in an unstable manner. 


17 Cracks in the shape of a disk lying in a plane are the only ones, aside from elliptical cracks, for which 
there is a known stress solution. 

18 A cylindrical rod of height, if, which undergoes a given small 5HH~ l when subjected to a uniaxial 
stress, will increase in radius by a 5R/R = v p 5HH~ l . Poisson’s ratio for an incompressible material is 0.5. 

19 The 7 , e for water is 77 ergs cm -2 . Most familiar solid materials have 7 se ’s well into the thousands. 
However, the 7 se for HMX crystals is only 50! 

20 If a solid body of height, if, is subjected to a small shear stress in an orthogonal direction to if by a 
pressure of P dynes cm -2 , and shears by 5L in that direction, then its shear modulus, a = PH{8L)~ l . 
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For the cases such as the end of data segments 20 and 22 (Figs. 4 and 5) which shows a 
very gradual onset of angular momentum transfer, it is necessary to show that crack growth 
can have a gradual, slow rise in velocity from zero, or a “creep” phase. Evans (1974) and 
Charles (1958) have observed that, when the stress is below the critical level of instability 
given above in Equation (1), cracks may grow at low speeds: 

» = »„(*/*,)* ( 2 ) 

where v is the crack growth speed, v R is the Rayleigh velocity, 0.93-y^/xp -1 , the speed of a 
surface (shear) wave in the material of density, p, K is the stress intensity factor at the crack 
tip, y'7rc(<7 — r), K\ is some stress intensity scale factor, and the optimum k is 12 for many 
brittle materials. The gradual behavior observed in data segments 20 and 22 (Fig. 5) may be 
the result of such a process. This formula holds only for a stress intensity factor, K, below 
a certain critical value, K 0 = y/2i tBjZ, where E is Young’s Modulus 21 , and %<, again, is the 
specific surface energy. For higher stress concentration factors, Freund (1990) has shown in 
theoretical studies that: 

v^v R (l-(Ko/K) 2 ), (3) 

and this equation may likely govern the crack speed behavior for all sudden glitches. A value 
for the scale of the stress concentration factor, K\, can be derived by setting the values and 
slopes of the v's for the two regions equal: 

Ki = K 0 ik/2)V k {l + 2/fc) 1 / 2+1 / fc . 

From equation 1, the critical radius for closed penny-shaped cracks is: 

_ 7T (2 - Up) JseP 
° cr 2(1 - up) (a -t) 2 ' 

and the corresponding fracture toughness is then: 

K 0 - - r) = 7T ^/(l — z//2)7 3e /^(l — u). 

And finally, the relative crack speed at which the two stress concentration factors are equal 
is: 

v/vr = 1/(1 + k/2), (7) 

which is 1/7 for k = 12. 


21 A body of length, JL, when subjected to a small tensile force per unit area, cr, will strain by an amount, 
5L — aLE “* 1 , where E is the Young’s Modulus (and 5L will drop to 0 when a = 0). The shear modulus, /z 
is related to E and Poisson’s ratio, up, by E = 2/z(l + up). 


(4) 

(5) 

( 6 ) 
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We can estimate 22 the critical crack radius from equation 5 by taking n = 2xl0 29 dynes 
cm -2 from Ruderman (1991b), the values for the displacement field near 0.05 cm from Baym 
and Pines (1971 - also as equations 1 & 2 in Link et al. 1998), a v P — 0.5 for incompressible 
matter, and amortizing this displacement over 1 km of crust to get a strain near 5xl0 -7 . 
Using the strain times a Young’s modulus of 3/jl = 6xl0 29 dynes cm -2 , yields a stress of 
3xl0 23 dynes cm -2 , or 3xl0 14 kBar, reassuringly close to yield stresses for normal, everyday 
metals scaled up to NS mid-crustal densities. For a realistic j se , we can scale up 10 4 ergs 
cm -2 , for a material with a mean density of 5, by lxlO 13 to match the density of the outer 
crust, and we can ignore, for the moment, the interfacial friction, r. Using equation 5 the 
resulting critical crack radius, for the onset of instability for these (hypothetical) conditions 
in the crust, is c cr = 5.5 cm. Of course this value will increase for any non-zero value for 
the interfacial friction, with a r of half of the 3xl0 14 kBar stress given above increasing 
^ to 22 cm. Thus the total growth needed per year for the glitches then approaches one 
meter, in agreement with the amount of pole migration required to account for the persistent 
increasing trend in \v\. 

The crust will also have some distribution of crack sizes, as occurs in the Earth’s crust, 
which is exponential over many orders of magnitude (Scholz 2002). In this case, there will 
be some probability that a crack of a certain size exceeding some limit will be present in 
the volume of interest. If the volume is not large enough, then there might well be no crack 
larger than that limit, and thus statistics enter into calculations of NS crust failure, just as 
it does in calculations of mechanical failure for other solid materials. 


22 Or calculate, as is done in Strohmayer et al. (1991) 
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Table 2. Master Pulse Profile (MPP) a ’ 6 


Constant Name 

Value 


a 

98,969. 

6 

b 

1,791,946. 


4> o 

0 . 

5050471985806 


0 . 

0104528014194 

Cl 

24. 

51180526856 

C 2 

295. 

1283580927 

c 3 

41. 

89721595187 

C 4 

1,168. 

592235964 

c 5 

111. 

6330356226 

c 6 

0 . 

3315320354509 

c 7 

3,145. 

824919010 


h MPP = aL(4> - (j> Q ) + 6; 0 . < (f> < 1.0 
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Table 5 , Log of Large Glitches from other Pulsars 
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Fig. 1. — The pulse profile for PSR J0537— 6910 from 311 observations taken during the 
RXTE campaign from Jan. 19, 1999, to Oct. 6, 2004. The parameters of the MPP curve 
drawn over the summed data are listed in Table 2. The vertical curve marks the midpoint 
of the pulse profile vs. height. 
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Fig. 2. — The frequency history for PSR J0537— 6910. The frequency measurements of 
the individual observations are shown with error bars. The (mostly horizontal) curves are 
the frequency histories derived from the much more accurate timing solutions. The right 
frame boundary labels the approximate cumulative gain of frequency across the glitches (see 
Fig. 8). Diamonds mark points with known phases from the timing solutions, but which 
were not included in them. Hexagons mark observations for which no timing solution exists. 
The horizontal dashed line is the baseline for the frequency residuals, and corresponds to 
62.037465 — 1.9922xlO~ 10 x(MJD — 51,400) x 86400 Hz. The oblique dashed line is an ap- 
proximation to the mean trend for the frequency residuals, and corresponds to 62.037545987 
- 1.9760x 10 -10 x (M JD - 51,400) x 86400 Hz. 




-46- 

(SJB3/C) ooz'pg am IB age guiunj, 


o 

00 

CO 


CO 

02 

02 

02 

a> 

. CD 

CD 

CD 







o’ o’ o’ o 


(puooas /z^ jshoubu) saiiba] 430 /Couanbsjji 

Fig. 3. — The v or spindown history for PSR J0537— 6910, derived from subsegments of 
at least three consecutive observations with a time span of at least 1,400,000 s (with the 
exception of the 17th data segment, which only spanned 14 days) , and at least one observation 
within 0.35 of the full time span of its center. The two dotted lines are spaced at 0.15 pHz 
s -1 and slope downward at 0.95xl0~ 21 Hz s -2 , or 0.15 pHz s -1 (five years) -1 , or a braking 
index, n, of —1.5. 
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Fig. 4. — The phase residuals to the timing solution fits for all 22 data segments, with the 
segment numbers listed just below the corresponding residuals. The solid points have not 
been included in the fits. The residuals for observations with no timing solutions have not 
been plotted. 
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Fig. 5. — (Lower) The timing residuals for data segment 7. The solid points have not been 
included in the fit. (Upper) The timing residuals for data segment 22. 
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Fig. 6. — The glitch amplitude vs. the time to the next glitch (see Table 4). The time bounds 
drawn for the points are equal (up and down) and each is a quarter of the sum of two time 
intervals bounding the data segment. The slope of the dashed line fitted to the points and 
through the origin is 6.444 days /uHz -1 , or 399.7 days ppm -1 . 
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Fig. 7. — The glitch amplitude vs. the following interval of stable timing behavior (see Table 
4). The slope of the dashed line fitted to the points and through the origin is 6.116 days 
/iHz -1 , or 379.3 days ppm -1 . The point for glitch 21 is very close to that for glitch 17 (24.9 
fiUz, 143.5 days, vs 24.3 /iHz, 143.1 days). 
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Fig. 8. — Cumulative time vs. cumulative glitch amplitude. The oblique line drawn from 
(0,0) to the bottom of glitch 21 has been offset from the first and last glitch times by ten 
days for clarity, and has a slope of 401.1 days ppm -1 or 6.468 days /zHz -1 . The top frame 
plots the predicted glitch onset time minus the actual glitch onset time (in days on the right 
hand scale, and converted to equivalent /j,Hz on the left hand scale). 








-54- 


uoi^^jSih SJBSiC 000‘0l/ 



(saAq) lua^q to ‘a^ay; 

Fig. 11. — (Lower) The magnetic pole migration rate vs. initial obliquity of PSR J0537— 6910 
for an age of 4,000 years, latency times (with no polar migration) from 2,000 to 3,600 years 
in steps of 400 years, and two values of the (intrinsic) braking index, n, 2.5 (solid) and 
2.0 (dotted/dashed). The four vertical lines of constant initial obliquity, or a 0 , intersect 
the curves at 22 migration rates, or d’s, for which time histories of v and v are plotted in 
Fig. 12 (see § 5.2.3). (Upper) Lines of constant present day magnetic obliquity, assuming 
a 4,000-year age for PSR J0537— 6910, for braking indices 2.5 (solid) and 2.0 (dashed) are 
plotted on the initial obliquity-time plane. See Fig. 12 and § 5.2.3. 
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Fig. 12.— Segments of 22 possible time histories of PSR J0537— 6910 for a present age of 
4,000 years, and for parameters from the intersections of ordinates of Fig. 11 representing 
initial magnetic obliquities of 0.1 - 0.4 radians in steps of 0.1 radians, with the lower migration 
rate segments of curves of constant pole-migration intervals from 2000 to 3600 years in steps 
of 400 years, and braking indices, ri, of 2.5 (solid) and 2.0 (dotted - see Fig. 11). The pulse 
frequency, v, histories are curves with values labeled on the left hand and bottom frame 
boundaries, first visible at the bottom center and continuing left/backward in time, while 
rising to the left-hand frame edge. The curves with the slope breaks are spindown histories 
and are labeled on the top and right-hand sides. 




